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SUMMARY 


Laminar  free  convection  from  a  non-lsothermal  horizontal  cylinder 
Is  analysed.  The  wall  surface  temperature  Is  assiuned  to  be  varied  In 
the  manner  of  transformations  are  devised 

and  employed  so  that  the  resulting  differential  equations  and  boundary 
conditions  are  free  of  the  parameters  a^  and  a^.  These  differential 
equations  ar4  solved  once  and  for  all)  solutions  to  the  original  equations 
for  any  particular  values  of  a^  and  a^  may  then  be  read  off  easily 
as  linear  combinations  of  the  numerical  solutions  given  here.  It  Is 
found  that  the  dependence  of  heat  transfer  from  a  horizontal  cylinder 
on  Prandtl  number  Is  practically  the  same  as  that  from  a  vertical  plate. 
Furthermore,  the  heat  transfer  is  greatly  Influenced  by  the  surface 
temperature  variations. 
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INTRCDUCTICN 

Laminar  free  convective  heat  transfer  from  a  cylinder  has  been 

stvLdled  theoretically  by  Hermann  [1]  and  by  Chlang  and  Kays  [2]. 

Both  references  report  the  solutions  of  the  pertinent  boundary  layer 

equations  for  the  case  of  an  Isothermal  vratll  surface,  the  first  for  a 

Prandtl  nximber  of  0.74  and  the  second  for  a  Prandtl  number  of  0.7. 

The  closely  related  problem  of  the  laminar  free  convective  heat  transfer 
f 

from  a  vertical  plate  has  been  analysed  by  Ostrach  [3]  for  an  Isothermal 
wall  and  with .a  wide  range  of  Prandtl  numbers.  When  either  the  plate 
surface  Is  non -Isothermal  or  the  surface  heat  flux  Is  prescribed,  the 
free  convection  problem  has  been  solved  approximately  by  Sparrow  [4] . 
Exact  formulations  and  solutions  have  been  reported  by  Sparrow  and 
Gregg  In  Refs.  [5]  and  [6]  for  uniform  heat  flux  and  non-isothermal 
surface  respectively.  A  search  of  the  literature  reveals  that  there  Is 
no  solution  available  for  the  prediction  of  free  convective  heat  trans¬ 
fer  to  or  from  a  cylinder  when  the  surface  Is  non -Isothermal.  As 
demonstrated  In  [6] ,  the  heat  transfer  from  a  non-lsothermal  plate  Is 
markedly  different  from  that  of  an  Isothermal  plate.  It  is  reasonable 
to  suppose  that  the  free  convective  heat  transfer  from  a  non-lsothermal 
cylinder  would  be  also  significantly  different  from  that  of  an  isothermal 
one.  Since  In  a  great  many  technical  applications,  the  cylindrical 
surface  from  which  heat  Is  being  transferred  Is  non-lsothermal.  It  Is 
the  purpose  of  this  paper  to  report  the  free  convective  heat  transfer 
from  a  non-lsothermal  cylinder.  The  surface  temperature  Is  assumed  to 
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take  the  following  form: 

(T^  -  V  =  -  Tj  [1  +  a^(f)2  +  (l) 

where  a^  and  a^  are  arbitrary  constants.  The  surface  temperature 
may  be  either  higher  or  lower  than  the  free  stream  temperature.  If  the 
surface  temperature  is  higher  than  the  free  stream  temperature,  x  must 
be  measured  from  the  lower  stagnation  point  (Figure  la) .  Oi  the  other 
hand,  if  the' surface  temperature  is  lower  than  the  free  stream  tempera¬ 
ture,  X  must  be  measured  from  the  upper  stagnation  point  (Figure  lb). 
When  the  coordinate  systems  are  taken  in  this  manner,  the  method  of 
analysis  and  the  results  for  the  heat  transfer  parameters  are  the  same 

for  either  T  >  T  or  T  <  ,  and  there  will  be  no  need  to  treat 

w  <»  w  <“ 

them  separately.  The  analysis  will  be  carried  out  for  the  case  of 
^  remembered  that  the  results  apply  to  both 

cases  depicted  in  Figure  1. 

GovER:;ii:a  difflreictial  equations 


The  equations  expressing  conservation  of  mass,  momentvun,  and  energy 
for  steady  laminar  flow  in  a  boundary  layer  on  a  horizontal  cylinder 
are  as  follows*: 


ar  9T  _  ji.  afr 

^ax  ^ay  Pr  gy2* 


(2) 

(3) 

(4) 


*  See  page  17  for  a  listing  of  the  definitions  of  the  various  variables 
Involved . . 
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The  density  has  been  considered  a  variable  only  .in  forming  the  buoyancy 
force  gp(T  -  T  )sln  Variations  of  all  other  properties  are  neglected. 
Viscous  dissipation  and  work  against  the  gravity  field  also  have  been 
omitted  from  the  energy  equation  (4)  since  they  are  negligibly  small. 

The  appropriate  boundary  conditions  aret 


y  =:  0  u=v=0  T  =  T^(x) 

y_«f  U->0  T-^T^ 


(5) 


Equation  (2)  can  be  automatically  satisfied  if  a  stream  function  is 
introduced  such  that 


(6) 


The  momentum  and  energy  equations  may  be  rewritten  in  the  following 
format 


M  M-  -  M  M-  =  M _ +  <)  sin  X 

f)  xti  fit)  X  nqri  ^ 


where 


-  X 

R 


=  [ - ^ ^  =  (Gr)^  i 


11 

M  =  [ 

0  = 


=  i 


wo 
T  -  T_ 


T  -  T_ 

wo  00 


(7) 

(8) 
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In  terms  of  the  new  variables,  the  bcnmdary  conditions  are 
fl  =  0  M  =  M-  =  M^=0,  <)=1+  a^x^  +  bl^x^ 

(10) 

^-*•00  M^-+0  V"*^0 


Similar  to  [l],  equations  (7)  and  (8)  together  with  the  bcundary 
conditions  (lO)  are  solved  by  a  perturbation  technique.  It  Is  assumed 
that  M  and  (()  take  the  following  format 


M  =  +  x^F^^Ctl)  +  5^F2(n)  +••• 

<)  =  G^(fl)  +  x^G^(ti)  +  x^G2(n)  +••• 


(11) 


When  equations  (ll)  are  substituted  Into  equations  (7)  and  (8) ,  and  equal 
powers  of  x  are  collected,  the  following  set  of  ordinary  differential 
equations  is  obtained; 


lit  II 

0  0  0 


1  " 
^  0 
Pr  o 


+  F  g'  =  0 
o  o 


(12) 


II  f  It  t  t  11 


5=0  =  » 


(13) 


^  “l  *  =  0 


in  II  1  f  11  11 

^2  ^0^2  -  ^^0^2  5^0  ^2  3^1^1  - 

.  .  ^2  -  ^1  ^  120  °o  =  ° 


|2 


(14) 


1  " 
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where  the  superscript  prime  (')  denotes  differentiation  with  respect 
to  ti  • 

Notice  that  in  obtaining  equations  (12)  and  (14),  sin  x  in 
equation  (7)  has  been  approximated  by  the  first  three  terms  of  its 
power  series.  For  ^  ^  ihe  errors  resulting  from  this 

truncation  are  about  0.5%  and  A%  respectively. 


The  boundary  conditions  are 


f|  -*•  «0 


F  =  F'  =  0  G  =  1 
o  .  o  . 


F '  -*•  0  ,  .  G  -*•  0 

o  o 

F£  0  Gl  •♦,  0 

F^,->0  Gg-^O 


(15) 


An  inspection  of  equations  (12)  to  (15)  reveals  that  the  present  problem 
involves  three  parameters,  Pr,  a^^,  and  a^.  For  each  combination 
of  Pr,  a^,  and  a^,  it  would  be  necessary  to  solve  three  pairs  of 
simultaneous  equations.  If  the  fluid  is  a  gas,  the  Prandtl  number  is 
relatively  constant.  However,  both  a^  and  a^  can  have  a  wide  range 
of  values,  and  hence  it  might  be  necessary  to  solve  equations  (12)  to 
(15)  for  a  great  many  combinations  of  a^  and  a^.  This  is  obviously 
impractical, if  not  impossible.  Consequently,  the  following  transforma¬ 
tions  (Eq.  16)  are  used  to  transform  the  original  equations  and  boundary 
conditions,  to  a  new  set  of  differential  equations  and  boundary  conditions 
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which  are  free  of  the  parameters  a^^  and 


^0  =  ^1 


^0=^1 


Fi  =  a^X^  +  X3 


(16) 


+  a^Xj  +  a^X^  +  =  a^Y^  +  a^Y^  +  a^Y^  +  Y^ 


The  new  differential  equations  together  with  their  boundary  conditions 
are 


f  2 

II I  *  II  t 

Xi  +  X^Xj^  =  X^  -  Y^ 

yJ  +  Pr  X^Y^  =  0 


X^(0)  =  X3^(0)  =  X^(»)  =  0 
Y^(0)  =  1  Y3^(«»)  =  0 


(I) 


ttf  tt  II  It 

^2  ^  "  "^1^2  *  ^^1^2  ■"  h 

X2(0)  =  X2(0)  =  X2(~)  =  0 
Y2(0)=1  Y2(«)  =  0 


(II) 


II I  It  II  n  1 

X3  +  Xj^X^  =  a^X^  -  3Xj^X3  -  Y3  + 

Y3  +  Pr[X3^Y3  -  2X^X3  +  =  0 

X3(  0)  =  X3(  0)  =  X3(«)  =  0 


(III) 


Y3(0)  =  Y3(o<.)  =  0 
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IM  II  I  I  It 

X4  +  =  6X^X4  -  5Xj^X4  - 

Y4  +  Pr[x^Y4  -  4X^Y4  +  5X4Y’]  =  0 


=  X4(0)  =  X4(«)  =  0 

Y4(o)  =  1  Y4(oo)  =  0 


lit  nit  ti  ,  •  «v 

Xj  +  X^Xj  =  6X^X5  -  5X^X5  +  3(X2  -  X2X2)  -  Yj 

Yj  +  PrtXj^Yj  -  AX^Y^  +  3X2Y2  -  2X2Y2  +  5XjY|]  =  0 
XjCO)  =  X^iO)  =  Xgloo)  =  0 
^5(0)  =  =  0 


(IV) 


\  +  =  6x^x^  -  5Xj^X^  +  3(  2X2X2  -  X2X2  -  X2X2)  - 

Xl  +  PrlXj^Y^  -  4xJ^Y^  +  3(X2Y2  +  X^Y^)  -  2(X2Y2  +  Y2X2)  +  5X^yJ[]  =  0 
X^(0)  =  X^(0)  =  X^(<»)  =  0  (VI) 

Y^(0)  =  Y^(<»)  =  0 


m  II  I  I  It  It  ^3 

^  +  X,x^  =  6X,X^  -  5X,X^  +  3(X3  -  X3X3)  -  ^ 

y!J  +  PrfX^Y^  -  4X^Y^  +  3X3Y3  -  2X3Y3  +  5X^Y^  =  0 
X^(0)  =  X^(0)  =  X^(~)  =  0 
Y^(0)  =  Y^(<»)  =  0 


(VII) 


Notice  that  equations  (l)  to  (VIl)  are  free  of  the  parameters  a^  and 
a2  and  hence  solutions  to  these  equations  are  applicable  for  all  values 
of  a^  and  a2. 
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Equations  (l)  to  (VIl)  have  been  solved  numerically  on  an  IBM  7090 
computer.  It  Is  noted  that  equations  (l)  are  non-llnear,  while  the  remainder 
of  the  equations  are  linear  (assuming  that  the  solutions  of  equations  (I), 
(ll),  and  (III)  are  now  known  one  after  the  other.  Thus  an  Integral  equa¬ 
tion  approach  was  thought  suitable  for  attempting  the  solution  of  equations 
( l) ,  while  Initial  value  methods  were  Indicated  for  the  remaining  equations . 


If  for  ease  In  writing,  one  defines  the  quantities 

n 

'  -5  X^(n)dri 

B^(tl)  =  e  ° 

B2(n)  =  [B^(»l)]^’' 


B3(r|) 


n  {[X'(fi)]^  -  Y  (n)| 
*  - 57^1 - 


00 

.(  B,(fl)BT(n).ir, 
0  ^ 

00 

■  5  B,(fi)dfi 
0 


a  system  of  Integral  equations  whose  solution  also  satisfies  (l)  and 
which  may  be  solved  by  a  method  of  successive  substitutions  Is 

n 

=  5  X'(n)dfi  (i7a) 

0 


5  B2('l)dTi 

=  h -  (IVb) 

5  B  (n)dn 

0 

n 

X4(n)  =  5  [B,(fl)  -  K]B^(fi)dt). 

0  ^ 


(17c) 
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Guesses  for  the  functional  values  X|(tl)  were  made.  Values  of 
and  were  then  obtained  by  (17a)  and  (l7b)  respectively.  These 
values  were  then  substituted  in  the  right  hand  side  of  (I7c)  so  as 
to  obtain  new  values  of  X|(ti).  (it  was  found  that  an  under-relaxation 
factor  of  approximately  0.3  would  insure  convergence  of  the  process 
when  reasonable  guesses  for  the  values  of  X|(ti)  had  been  given  originally.) 

After  equations  (l)  had  been  solved  fairly  accurately  by  the  method 
just  outlined,  one  had  in  particular  accurate  approximations  for  X^(0) 
and  Y^(0).  Then  even  on  the  non-linear  equations  (l),  initial  value 
methods  became  very  attractive.  Actually,  the  Runge-Kutta  method  was 
used  to  obtain  final  results  which  aie  better  than  those  obtained  by 
the  integral  equation  approach}  they  are  better  mainly  because  it  was 
feasible  to  use  a  much  smaller  value  of  Aii  in  the  initial  value  method 
than  in  the  integral  equation  method.  Finally  the  Runge-Kutta  method 
was  used  also  to  solve  equations  (ll)  through  (VIl) . 

For  the  problems  in  which  we  have  been  interested,  the  Prandtl 
munber  generally  ranges  from  0.7  to  1.0.  In  Tables  I  through  VII, 
tabular  results  are  given  for  the  cases  Pr  =  0.7  and  Pr  =  1.0.  In 
accordance  with  equations  (16),  linear  combinations  of  the  tabular 
results  will  be  solutions  of  the  original  differential  equations  (12)- 
(14) }  the  particular  linear  combination  desired  will  depend  on  the 
values  of  a^  and  a.^  in  the  boundary  conditions  (15) . 
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VELOCITY,  TEMPERATURE  AND  HEAT  .TRANSFER 


The  velocity  and  temperature  distribution  in  the  boundary  layer 
can  be  calculated  by  the  equations 

=  xX^  +  +  X^]  +  +  a^X^  +  +  X^] 

T  -  T 

^ ^  =  Y^  +  x2[a^Y2  +  Y3]  +  ^^[agY,  +  a^Y^  +  a^Y^  +  Y^] . 

wo  00 

f 

The  dimensionless  heat  transfer,  the  Nusselt  number,  can  bo  computed 
from  equation  (20) . 


uR 

Gr^v 


(18) 

(19) 


i!!^=  JlS_=  _  {y.(o)  +  x^fa  Y'(0)  +  Y'(0)] 

Gr*  kGr* 

+  x^[a2Y|(0)  +  a^Y^(O)  +  a^Y^Co)  +  Y^(0)] 


(20) 


where 


For  Pr  =  0.7, 

=  .37023  +  [.75688aT  -  .01609] (|)^ 

Gr^  ^  ^ 

+  [.92847a2  +  .09471a^  -  .02885aj^  -  . 00009] (|)^.  (2l) 

For  Pr  =  1.0, 

^  =  .42143  +  [.85200a.  -  .01861](|)^ 

Gr^  ^  ” 

+  [l.0411a2  +  .10702a^  -  .03305aj^  -  .0001l](|)^.  (22) 


DISCUSSION  OF  RESULTS 


Velocity  and  Temperature  Profiles  The  velocity  and  temperature 

profiles  as  computed  from  equations  (18)  and  (19)  for  Pr  =  .7  at 

^  =  1  are  plotted  in  Figures  2  and  3*  It  is  evident  from  Figures  2 
K 

and  3  that  an  increase  in  wall  surface  temperature  will  Increase  the 
maximum  velocity  in  the  boundary  layer  and  will  increase  the  absolute 
value  of  the  temperature  slope  near  the  wall, 
f 

For  a^  =  a^  =  0  (isothermal  case)  and  Pr  =  0.7,  a  comparison 
of  [2]  and  the  present  calculation  reveals  that  the  results  listed  in 
Tables  I  and  III  agree  with  those  in  [l]  very  well.  However,  there  are 
significant  discrepancies  between  the  results  of  Table  VII  and  those  in 
[2].  When  the  initial  values  as  given  in  [2]  were  used  in  the  integra¬ 
tion  of  equations  (VIl),  it  was  found  that  the  conditions  at  "free 
stream"  were  not  satisfied.  Consequently,  it  is  believed  that  the 
values  listed  in  Table  III  of  [2]  (corresponding  to  Table  VII  of  the 
present  paper) are  in  error.  However,  since  the  absolute  values  are 
quite  small,  this  error  may  not  have  a  significant  effect  on  the  overall 
results . 

A  comparison  of  velocity  and  temperature  profiles  in  the  boundary 
layer  was  made  in  [2]  between  the  results  of  [l],  [2],  and  the  existing 
experimental  data  for  ^  “  J  I'*'  concluded  that  they  agree 

very  well  with  one  another.  This  same  conclusion  is  directly  applicable 
to  the  present  case  of  a^  =  a2  =  0.  As  stated  before,  we  know  of  no 
previous  theoretical  or  experimental  results  for  the  case  of  a  non- 
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Heat  Transfer  Figure  4  shows  the  effect  of  the  Prandtl  nvunber  on 
isothermal  wall  heat  transfer.  When  the  Prandtl  number  is  increased 
from  0.7  to  1,  the  heat  transfer  is  increased  by  about  13  to  14 
percent,  depending  on  the  location.  It  is  interesting  to  compare  the 
effect  of  the  Prandtl  number  on  heat  transfer  from  a  vertical  plate  to 
that  from  a  cylinder.  It  was  derived  by  Eckert  [7]  and  confirmed  by 
Ostrach  [3]  that  the  effect  of  the  Prandtl  ntunber  on  free  convective 
heat  transfep  from  an  Isothermal  vertical  plate  may  be  approximately 
computed  from 

NU  NU  ,  nrn  1  i 

^  ^  .952  +  Pr^ 


If  the  Prandtl  number  is  Increased  from  0.7  to  i,  the  heat 
transfer  as  computed  from  equation  (23)  would  be  increased  by  about 
13  percent.  In  view  of  this  close  agreement  on  numerical  results  and 
in  view  of  the  fact  that  the  heat  transfer  from  a  vertical  plate  and 
a  horizontal  cylinder  are  intimately  related  [7] ,  it  is  recommended 
that  in  the  absence  of  heat  transfer  results  for  the  Prandtl  numbers 
different  from  0.7  and  1,  equation  (23)  be  used  to  compute  the  free 
convection  heat  transfer  for  different  Prandtl  numbers. 

The  effect  of  a  non-isothermal  wall  on  heat  transfer  may  be 
computed  from  the  following  equation; 


r  q  ll 

T 

H 

iso 

/T  -  T  \ 
wo  W  1 

L(T„  -  TjGrtJi 

-  TjGrl 

'  W 

(24) 
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where  the  subscript  iso  denotes  a  quantity  corresponding  to  an 

^"'^00 

isothermal  wall  (a^  =  a^  =  O) .  The  quotient  {— - is  prescribed, 

while  the  remaining  quantities  on  the  right  hand  side  of  equation  (24) 
can  be  readily  calculated  by  use  of  equations  (2l)  and  (22) .  The 
results  of  calculation  for  Pr  =  0.7  and  two  combinations  of  a^  and 
a^  are  shown  in  Figure  5-  It  is  seen  from  Figure  5  that  the  wall 
temperature  variation  has  a  significant  effect  on  heat  transfer.  For 
a^  =  Sg  =  1  fat  the  location  ^  ~  2  (about  115°) ,  the  heat  transfer 
is  triple  the  value  of  that  corresponding  to  the  isothermal  cylinder. 
Clearly,  this  indicates  that  the  heat  transfer  in  free  convection  from 
or  to  a  non-lsothermal  cylinder  may  not  be  reliably  calculated  from 
solutions  corresponding  to  an  .isothermal  wall.  This  conclusion  was 
drawn  by  Sparrow  and  Gregg  [5]  for  the  case  of  a  vertical  plate. 

Comparison  of  Isothermal  Heat  Transfer  As  stated  previously,  the  free 
convection  heat  transfer  from  an  isothermal  horizontal  cylinder  was 
first  solved  by  Hermann  [l].  It  is  therefore  advisable  to  compare  the 
present  heat  transfer  results  with  those  of  [l] .  Since  Hermann's 
results  are  valid  for  Pr  =  0.74,  it  is  necessary  to  modify  his  results 
so  that  the  comparison  may  be  made  based  on  the  same  Prandtl  number 
(l.e.,  Pr  =  0.7).  Using  equation  (23),  it  is  found  that  the  correction 
factor  is  0.976,  l.e., 

(-^)  =  0.976(^) 

Gr^  Pr=0.7  Pr=.74 


(25) 
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Using  equation  (25) )  the  heat  transfer  results  as  given  by  Hermann 
and  those  of  the  present  calculations  for  Pr  =  0.7  are  listed  as 
follows! 


Q  I 

0 

TT 

? 

TT 

3 

IT 

2 

2iT 

3 

5tt 

T 

fl, 

Nu 

(RefQ) 

0.376 

0.372 

0.356 

0.329 

0.288 

0.227 

0.178 

0 

O) 

Present 

0.370 

0.366 

0.352 

0.330 

0.298 

0.256 

— 

— 

Ratio  ^ 

.984 

.982 

.992 

1.004 

1.035 

1.127 

B 

B 

It  is  seen  from  the  above  table  that  the  heat  transfer  results 
agree  with  each  other  to  within  3*5  percent  for  ^  Near  x  =  0, 

the  present  results  are  exact,  and  hence  Hermann's  calculations,  after 
being  corrected  for  the  Prandtl  number  effect,  underestimate  the  heat 

X 

transfer  by  about  2  percent.  For  5-  >  -5-  ,  Hermann's  results  are  more 
accurate  since  the  perturbation  solution  used  in  this  paper  breaks  down 
for  —  approaching  tt.  In  order  to  use  the  present  solutions  for  heat 
transfer  calculation  at  ^  it  is  suggested  that  the  following 

formula  be  used: 


5>2it/3 


=  Y( 


T  -  T 
w  ® 


x=2tt/3 


(26) 
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The  correction  factor  Y  can  be  derived  from  [l]  and  its  values 
are  as  follows: 

IT 

0 

The  correction  factor  Y  listed  above  is  strictly  for  an  iso¬ 
thermal  cylinder.  From  physical  considerations,  the  boundary  layer 
is  very  thick  near  x  =  tt  and  hence  heat  transfer  is  very  small  even 
when  the  wall  is  non -isothermal.  Therefore,  before  more  reliable 
results  are  available,  it  is  suggested  that  the  above  correction 
factor  be  also  used  for  the  non-isothermal  case.  Hence  one  can  first 
compute  the  heat  transfer  for  x  =  -j-  by  use  of  equations  (ll)  and  (12) 
and  then  obtain  the  heat  transfer  for  x  >  -j-  from  equation  (26)  to¬ 
gether  with  the  correction  factor  Y  given  in  the  table. 

In  many  practical  problems  the  flow  in  the  boundary  layer  far 
downstream  of  the  stagnation  point  may  be  turbulent.  In  this  case, 
the  local  heat  transfer  must  be  estimated  from  an  appropriate  turbulent 
heat  transfer  formula. 


X 

R 


2tt 

3 


5" 

T 


11 

12 


.788 


.620 


CONCLUSIONS 

The  boundary  layer  differential  equations  for  free  convective  flows 

over  a  horizontal  cylinder  are  solved  by  a  perturbation  method.  The 

wall  surface  temperature  is  assumed  to  be  varied  in  the  manner  of 

aT(^)^  +  ao(“)^»  A  special  transformation  is  used  so  that  the  solutions 
1  R  2  R 
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are  Independent  of  the  parameters  a^^  and  a^,  and  hence  results  are 
applicable  for  all  values  of  a^^  and  a^.  It  Is  found  that  the  functional 
relationship  between  heat  transfer  and  the  Prandtl  number  for  a  cylinder 
Is  essentially  the  seme  as  that  for  a  vertical  plate.  JVirthermpre,  the 
wall  temperature  variation  has  a  significant  effect  on  heat  transfer  and 
the  heat  transfer  from  a  non -isothermal  wall  m''y  not  be  computed  from 
solutions  corresponding  to  an  Isothermal  wall. 
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^l’^2 

c 

P 

F,  G 
Gr 

g 

h 

k 

M 

Nu 

Pr 


q 

R 

T 

T 

1 

T 


wo 

00 

u,  V 

X,  Y 
X,  y 

X 


p 

n 

p 


NOMENCLATURE 

Wall  Surface  Temperature  Parameters,  Eq.  (l) 

Specific  Heat  at  Constant  Pressure 
Eqs.  (11) 

Grashof  Number,  gp(T^^  -  T^R^/v^ 

Gravitational  acceleration 

Heat  Transfer  Coefficient,  h  =  q/(T^^  -  T^ 

Thermal  Conductivity 

Dimensionless  Stream  Rinction,  Eq.  (9) 

Nusselt  Number,  [q/(T^^  -  T^](R/k) 

Prandtl  Number,  pc ^/k 
Surface  Heat  Flux,  -  (k  5T/8x)^ 

Radius  of  Cylinder 
Temperature 
Surface  Temperature 

Surface  Temperature  at  Lower  Stagnation  Point 
Ambient  Temperature 

Velocity  Components  Along  x  and  y  Directions,  Respectively 
Dependent  Variables,  Eq.  (l6) 

Coordinates  (Fig.  l) 

Dimensionless  Coordinate,  Eq.  (9) 

Expansion  Coefficient 
Dimensionless  Coordinate,  Eq.  (9) 

Dynamic  Viscosity 
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V  Kinematic  Viscosity,  (i/p 

p  Density 

(|)  Dimensionless  Temperature,  Eq.  (9) 

'll  Stream  Rmction,  Eq.  (6) 

Subscripts 

“  at  infinity 

iso  for  isothermal  wall 

0  at  X  =  0 

w  at  wall 

wo  wall  surface  at  x  =  0 

Superscript 

'  differentiation  with  respect  to  q 
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Fig.  1.  Coordinate  Systems 


?0 


n 


Fig  2.  Velocity  Profiles 


Fig  3.  Temperature  Profiles 
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R 

Fig  4.  Effect  of  Prandtl  Number  on  Heat  Transfer 


J 

Fig  3.  Effect  of  Non  isothermal  Wail  on  Heat  Transfer 


72 


TABI£  I  SmjffirNS  OP  EWATIONS  (I) 


Pr  ■  0,7 


Pr  • 

1.0 

0 

1  *'l 

mm 

1  ’f'l 

.859 

54 

1.0000  1 

-.37023 

.81 Toil 1.0000 


.62622 

.9157I) 

.53flli6 

.87368 

.1)5626 

.03176 

.37902 

.79007 

.  3092*1 

.71)073 

.21)1)56 

.707O5 

.39301  .13250  .620o1i 


1 

I  .32705 

.111803 

.06597 

.60286  -.33665 

.30378 

.Ol270 

.55172 

2 

.37219 

.11257 

.51518 

SB 


.59555 


.63709  .'11755 


.67901  i  .liOitTO  -.13'<'*.7 


.50721  I  .39100  -.097001  .30233 


.5'*500  .30126  -.III1O7 


.36910  -.12703  .32lt06 


.35595 


2000 


.9O167  .205't0 


-.11192 

.15650 

-.13110 

.13217 

-.12301 

.11172 

-.09960 

.O7OO5 

.92120 

.19120 

.95755  1 

.I6O95 

.10731  .0053l<  .Oltli75 


1.3301t  .OO95O 


i.sI'I'tI  .00300 


1.3*173 


1.3‘i00 


1.3'> 


1.3'tO 


1.1200 


1.1559 


.(10160 -.00151  1.1595  1.00212 


U.608 


.OOII6I-.OOI35 


ifv  |%o  1 Ico  I ov  I  |o  I  l>o 


TABLR  tt  nOLWIOHn  Of  KQIIATIOfia  (II) 


Pr  •  0.7 


Pr  ■  1.0 


.ii6n>«8  1.0000 


.0031S  .0)1198 


.00809  .  07)179 


.01686  .  09936  I  .?0760  I  .777591 -.71?72 


.02778  .1166)1  .13959 


.03998  I  .12760  I  .08113  |  .6)il01 


.00205 

.039)19 

.00760 

.06991 

1.0000 

-.85800 

.91506 

-.8)i)ik0 

.83155 

-.82380 

[B!TnBB!3BfSBIi^!3BSS3WyHT!!| 


-.61356 

.0)1896 

.12070 

.0181 

)i  .53108 

-.57503 

.06106 

.12056 

-.0196 

8  .)i67l)i 

>.56)176 


.08)138  I  .110861  -.072581  .35)i18|-.51510 


.11836 


.10893 


.0983)1 


.08705 


.1)1)108  .075)10 


.28006 


.2)1379 


.09507  .10272 


.10)187 


.11367 


.121)10 


.12801 


.13352  .0)1961 


.13795  .03896 


.1)1133  .02881 


.1)1373  .01930 


.1)1522  .01052 


-.  >16652 


-.I1I967 


-.37501 


-.33290 


-.29356 


-.2571k 


-.22370 


-.1932)1 


-.16569 


-.lk097 


.11036 


5.0  1.079)1)1 


.06117 


iSBBBtl 


.13020 


.12310 


-.00735  -.01837  •11568 


-.01021  -.01063  *10807 


.00)152  -.01175  -.00508  . 10066 


.008)10  I-.01235 


BBBilBBBlBPCPlI 


-.018961  -.011511 -.01888 


-.0lkl9 


-.oom  -.015161 -.00167 


.00396^  -.01503 


.00809  -.0l)i21|  .00525 


.010861  -.01301 


.01105  I -.012321 -.00139 


.  00303".  0)1767 


.00)178 1 -.00215 .00163 


.00237  I -.000991  .00078 


00109 


.00050  -.00019  .00016  .036)17 


.00009  I -.0000)1  .000031.03621 


.00355  -.000831  .00083 


.00160  -.00030  .00031 


.00068  -.00011  .00011 


.00028  -.0000)1 


0 


TAWi  III  SOLWriOHS  Of  EWTIOftO  (III) 


Pr  ■  O.Y  Pr  ■  1.0 


3  X-3  X-j 


-.OlItTlI-. 05033 


-.035l*o| 

1  .0061*0 1 

1  .01570 

-.02367 

-.03219 

.00738  1 

.OlOOl 

iksssibsbihbtswshi 


1.0 

[M«aKBa 

.02279 

1 .0201*5 

1  .00198 

-.01931  .0225? 


-.01703  .03202 


-.011*73  .02300 


-.012l»li  .02250 


-.01022  .02101 


-.OO0O9  .02070 


.01967  .00723 


.02031  .00551 


.02077  .OO3OI 


.02107  .00215 


.02121  .00057 


.OlOl**  .02035 


.0151*6  .01973 


.01003 


.00530  I  .007191  .0l6l*8|-.006l*7 


-.00671 


-.00620 


.00730 

.00306 

.01381 

.00775 

.00151 

.0121*0 

-.00607 

.01955 

-.0021*1* 

.01671 

.00060 

.01366 

.00303 

.01065 

.001*87 

.00707 

.00620 

.0051*1 

.00706 

.00331 

.00751*1 

.00159 

.00772 

.00022 

.00766 

-.OOOO2 

.00591*  -.00253  .00521 


.00355  -.00200  .00251 


00127 


.00093  -.00060!  .00052 


.0001*1*  -.00031*  .00023 


.00020  I -.00016  .  00010  I  .00008 


.0071*1 

-.00150 

.00676 

.00502 

-.00263 

.00262 

.00273 

-.00185 

.00109 

.00132 

.0001*0 

.00025 1 -.00022  I  .00005 


0 


TABt/.  IV  BOLOnOWS  OF  TOUATIOHS  (IV) 


Pr  -  0.7 


Pr  ■  1.0 


.Ol*?15 


.05208  .09802  -.0352I< 


.08168  .09323 


.07086 


.07889 


.08825  .08893  |-.09l<91 


. 10258 


.10807 


.10867  .  02163 


.110>il  .013>>0 


.11137  .00585 


.00098 


1.0000 


.0017<i  .03325  .28815  .89839 


.00635  I  .05782 1  .20328  I  .7955 


.01301  I  .07lt.30l  .132ltl  .89952 


.02100  .081(1)9  .0733>)  .60968  -.881(33 


.02973  .089311  .02533 


.03872  .08990 


.01(780  .08711 


.08182 


.06390  I  .07k7‘i 


.07097  .0881(8 


.07717  .0571(9 


.0211)2 


.091)58  .013>il( 


.0081? 


.38359  -.8kl5<) 


|1(802 

.0821(6 

.01(822 

iSSSlI 

.11)872 

10506 

.08882 

.03897 

.11853 

,03588 


.023V5 


.01333 


.00213  -.011)1)1 


-.011(63 


.01056  -.01333  .001(93 


IS3E9I 


IBBTlHi 


.021(01 


.01530 


.00116  I  .01136 


iS339i 


.00862  I-.OO397  .00382 


.001)61)  -.0011(8  .00150 


.00217  I-.OOO51I  .00051) 


iBcinissisi 


16.0  1  .01070 


.00902 


.00675  I-.00016  1 .00015 


00002  I  .00880  .00003  .00002 


.00857 


.00855 


.000081  .00006 


.00002 


TAHU:  V  ;!n|.imntlfi  ok  K^llATIOIin  (V) 


?7 


TABLK  VI  SOLUTIOHS  OF  EQDATIOHC  (VI) 


TAJiii;  VII  noLin'ion;;  or  kwatiohs  (vii) 


10.0 

12.0 

lli.O 


R]-;KKRf:N(:h;s 
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